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Based on a sample of 2.25 x 10 8 J/ip events collected with the BESIII detector at BEPCII, a full 
partial wave analysis on J/ip — i- 77777 was performed using the relativistic covariant tensor amplitude 
method. The results show that the dominant ++ and 2 ++ components are from the /o(1710), 
/o(2100), / (1500), /2(1525), / 2 (1810) and / 2 (2340). The resonance parameters and branching 
fractions are also presented. 

PACS numbers: 13.25.Gv, 14.40.Be, 13.40.Hq 

I. INTRODUCTION 

Our present understanding of the strong interaction is based on a non-Abelian gauge field theory, Quantum Chro- 
modynamics (QCD), which describes the interactions of quarks and gluons; it also predicts the existence of new types 
of hadrons with explicit gluonic degrees of freedom (eg. glueballs, hybrids and multi-quarks) [1-5]. These uncon- 
ventional states, if they exist, will enrich the meson spectroscopy greatly and shed light on the dynamics of QCD. 
According to lattice QCD predictions [6, 7], the lowest mass glueball with J PC = ++ is in the mass region from 
1.5 to 1.7 GcV/c 2 . However, the mixing of the pure glueball with nearby qq nonet mesons makes the identification 
of the glueballs difficult in both experiment and theory. Radiative J/ip decay is a gluon-rich process and has long 
been regarded as one of the most promising hunting grounds for glueballs. In particular, for a J/ip radiative decay 
to two pseudoscalar mesons, it offers a very clean laboratory to search for scalar and tensor glueballs because only 
intermediate states with J PC = even ++ are possible. An early study of J/ip — > 77777 was made by the Crystal Ball 
Collaboration [8] with the first observation of /o(1710), but the study suffered from low statistics. 

In this paper, the results of partial wave analysis (PWA) on J/ip —> "/r)T) are presented based on a sample of 
2.25 x 10 s J/ip events [9] collected with the Beijing Spectrometer (BESIII) located at the upgraded Beijing Electron 
and Positron Collider (BEPCII) [10]. 

II. DETECTOR AND MONTE CARLO SIMULATION 

The BESIII detector, described in detail in Ref. [11], has an effective geometrical acceptance of 93% of 47r. It contains 
a small cell helium-based main drift chamber (MDC) which provides momentum measurements of charged particles; a 
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time-of-flight system (TOF) based on plastic scintillator which helps to identify charged particles; an electromagnetic 
calorimeter (EMC) made of Csl (Tl) crystals which is used to measure the energies of photons and provide trigger 
signals; and a muon system (MUC) made of Resistive Plate Chambers (RPC). The momentum resolution of charged 
particles is 0.5% at 1 GeV/c in a 1 Tesla magnetic field. The energy loss (dE/dx) measurement provided by the MDC 
has a resolution better than 6% for electrons from Bhabha scattering. The photon energy resolution can reach 2.5% 
(5%) at 1.0 GeV in the barrel (endcaps) of the EMC. And the time resolution of TOF is 80 ps in the barrel and 110 ps 
in the endcaps. 

Monte Carlo (MC) simulated events are used to determine the detection efficiency, optimize the selection criteria, 
and study the possible backgrounds. The simulation of the BESIII detector, where the interactions of the particles 
with the detector material are simulated, is GEANT4 [12] based. The J ftp resonance is produced with KKMC [13, 14], 
while the subsequent decays are generated with EvtGen [15]. The study of the background is based on a MC sample 
of 2.25 x 10 8 J/ijj inclusive decays which are generated with known branching fractions taken from the Particle Data 
Group (PDG) [16], or with lundcharm [17] for the unmeasured decays. 

III. EVENT SELECTION 

In this analysis, the 77 meson is detected in its 77 decay. Each candidate event is required to have five or six good 
photons and no charged tracks. The photon candidates are selected from the showers in the EMC with deposited 
energy in the EMC barrel region (| cosd\ < 0.8) and EMC endcap region (0.86 < | cos#| < 0.92) greater than 25 MeV 
and 50 MeV, respectively. The energy deposit in nearby TOF counters is included to improve the reconstruction 
efficiency and energy resolution. 

To suppress the background events with ir° (eg. J/ip — > 77r°7r°), the events that satisfy |M 77 — m w o| < 0.015 GeV/c 2 
are removed, where M 77 is the invariant mass of any pair of photon candidates and is the nominal ir° mass [16]. 
Then a four-constraint kinematic fit (4C), imposing energy-momentum conservation, is performed under the J/iJj —> 57 
hypothesis to reduce background events and improve the mass resolution, and xlc * s required to be less than 50. If 
the number of selected photons is larger than five, the fit is repeated using all permutation of the photons, and the 
combination with the smallest \\c * s selected. 

To distinguish the photons from 7/ decays, a variable 5, defined as S = \/{M lll2 — m v ) 2 + (A/ 7374 — to?,) 2 , is 
introduced, and the combination with the minimum value of S is chosen. The scatter plot of the invariant mass of 



one 77 candidate versus the other is shown in Fig. 1 (a), where the decay J/ip — » 77777 is clear. In order to select a 
clean sample, both M 7l72 and M 7374 are required to be in the 77 mass region, |M 7l72 (M 7374 ) — m v \ < 0.04 GeV/c 2 , 
with m v the nominal 77 mass [16]. The mass resolution for m v is about 10 MeV/c 2 . 

MC study shows that after the above selection, about 5.3% of events have a mis-combination of photons, which 
mainly occurs between the radiative photon and one photon from an 77. Therefore, candidate event must have only one 
combination with 6 < 0.05 GeV/c 2 to remove these events, which reduces the fraction of events with a mis-combination 
of photons to be 0.8%. 

After that, clear diagonal bands, which correspond to the structures observed in the 7777 invariant mass spectrum, 
can be seen in the Dalitz plot for the selected J/x/j — > 77777 candidate events (Fig. 1 (b)). A further requirement on 
Mj V , \M 1V — to^I >30 MeV/c 2 , is used to reject background events from J/tp —> <fii](4> — > 777). Fig. 1 (c) shows the 7777 
invariant mass spectrum of the surviving 5460 events after the event selection. 
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FIG. 1: (a) The scatter plot of M 7172 versus M 7374 after selecting candidates with the minimum S (two entries/event). The 
two-dimensional 77 sidebands are framed in regions A and B. (b) Dalitz plot (two entries/event), before <j> 777 veto, (c) The 
invariant mass spectrum of 7777. The dots with error bars are data, and the shaded histogram is background estimated from 77 
sidebands. 

Potential background events have been studied using both MC and data samples. The non-77 events are determined 
from the two dimensional 77 mass-sidebands, 0.07 GeV/c 2 < |M 77 — m v \ < 0.15 GeV/c 2 , which are defined as frame A 
and B in Fig. 1 (a). The shaded histogram in Fig. 1 (c) shows the sideband events normalized according to J/tp —> 57 
phase space MC. In addition, the background events are studied with a MC sample of 225 million J/ip inclusive events, 
and the main background events are found to be from J/ip — > 7777r°7r° and 77^77. In this paper, the background events 
estimated from 77 mass-sidebands, corresponding to a background level of 6%, are used in the partial wave analysis 
below. 
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IV. PARTIAL WAVE ANALYSIS 
A. Analysis method 

With GPUPWA, a Partial Wave Analysis framework harnessing GPU parallel computing [18], a PWA was performed 
to disentangle the structures present in J/ip — > 7?/?/ decays. The quasi two-body decay amplitudes in the sequential 
decay process J/tp —> "yX,X — > rjrj are constructed using covariant tensor amplitudes described in Ref. [19]. For J/ip 
radiative decay to mesons, the general form for covariant tensor amplitude is: 

A = ^(mi)e*(m 2 )^ = ^(m 1 ) e :(m 2 )I] l A i rj^, (1) 

where ip^imi) is the J/ip polarization four- vector, e v {rn-i) is the polarization vector of the photon and XJf v is the 
partial wave amplitude with coupling strength determined by a complex parameter A$. The partial wave amplitudes 
Ui for the intermediate states used in the analysis are constructed with the four-momenta of daughter particles, and 
their specific expressions are given in Ref. [19]. 

For an intermediate resonance, the corresponding Breit-Wigner propagator is described by a function: 

BW ^ = 1 —T7F> ( 2 ) 

M z — s — iMl 

where s is the invariant mass-squared of daughter particles, and M and T are the mass and width of the intermediate 
resonance. 

The relative magnitudes and phases of the amplitudes are determined by an unbinned maximum likelihood fit. The 
resonance parameters are optimized by the scan method: repeating the fits iteratively with various masses and widths 
until the optimized likelihood values converge. For the tensor states, the relative phases between three amplitudes 
for a certain resonance are theoretically expected to be very small [20] ; therefore the relative phases are set to zero in 
the fit so as to constrain the intensities further. 

The basis of the likelihood fitting is that a hypothesized probability density function (PDF) would produce the 
data set under consideration. The probability to observe the event characterized by the measurement £ is : 

where e(£) is the detection efficiency and w(£) = ^ is the differential cross section, and d(f> is the standard clement 
of phase space. The full differential cross section is: 

(1(7 



^^| 2 = \A{Vn) + A(0 ++ ) + A(2 ++ ) + A(4 ++ ) + ...| 2 , (4) 



where A(J PC ) is the full amplitude for all resonances whose spin-parity are J PC , and A(Vrj) is the contribution of 

the sequential decay processes such as J/ip —> (pi] — > 77777. / d£cj(£)e(£) = a' is the measured total cross section. 

The joint probability density for observing the N events in the data sample is: 

N N (*z.\. e (p.\ 

£ = U P &) = Ii-^ M - (6) 

i=l i=l a 

For the technical reasons, rather than maximizing C, iS=-ln£ is minimized, i.e., 

i=l i=l 

for a given data set, the second term is a constant and has no impact on the determination of the parameters of the 
amplitudes or on the relative changes of S values. So, for the fitting, ln£ defined as: 

N ,da_\ N , 

m£=x;in(^)=E i »©'- JVW ' (7) 

i=l i=l 

is used. The free parameters are optimized by FUMILI [21]. The measured total cross section a' is evaluated using 
MC techniques. Namely, a MC sample of N gen is generated with signal events that are distributed uniformly in 
phase space. These events are subjected to our selection criteria and yield a sample of iV acc accepted events. The 
normalization integral is computed as: 

r 1 Nacc a 

The background contribution is estimated with 77 sidebands. In the log likelihood calculation, the likelihood value 
of 77 sidebands events are given negative weights, and are removed from data since the log likelihood value of data is 
the sum of the log likelihood values of signal and background events, i.e., 

S = -(ln£rH UH - In Leg). (9) 

The number of the fitted events Nx for an intermediate resonance X, which has Nyy independent partial wave 
amplitudes A;, is defined as: 

N X = ™ ■ N', (10) 
a 

where N is the number of selected events after background subtraction, and 

N acc N 



face 1 *w 



ax ~ N 

Jv acc , 

k 3 



is the measured cross section of the resonance X and is calculated with the same MC sample as the measured total 
cross section ct'. 

The branching ratio of J/ip — > jX, X — > rfC) is calculated with: 

B{J/i> ->• lX 77777) = — , (12) 

Nj/tp ■ e x ■ o,_> 77 

where the detection efficiency ex is obtained by the partial wave amplitude weighted MC sample, 

_ ox Ef cc IEf w (^) fc | 2 m , 
£x °F E^E^OW' 1 j 

The statistical errors for masses, widths and branching ratios in a PWA are defined as one standard deviation from 
the optimized results, which corresponds to a change, 0.5, of the log likelihood value for a specific parameter. In this 
analysis, the changes of log likelihood value and the number of free parameters in the fit with or without a resonance 
are used to evaluate the statistical significance of this resonance. 



B. PWA results 



In this analysis, all possible combinations of ++ , 2 ++ , 4 ++ resonances listed in the PDG summary table [22] are 
evaluated, and the fitted components with statistical significance larger than 5. Oct are kept as the basic solution. The 
contribution from 4 ++ (/ 4 (2050)) with a statistical significance of 0.4ct is ignored. There are six resonances, /o(1500), 
/o(1710), / (2100), /a(1525), / 2 (1810), / 2 (2340), as well as 0++ phase space and J/ip -> cjyq included in the basic 
solution. Although most of the J/ip — > <pr) events have been rejected by the above <p mass window requirement, 
J/ip — > <p>r\ is included in the PWA to evaluate its impact from the interference between the tail of <p an d other 
components from J/ip — > ^X(X — > 7777). The masses and widths of the resonances, branching ratios of J/tp radiative 
decaying to X and the statistical significances are summarized in Table I. 

The comparisons of the r\r\ invariant mass spectrum, cosO^, cos(9 7 and (p^ distributions between the data and the 
PWA fit projections (weighted by MC efficiencies) are displayed in Fig. 2 (a), (b), (c), and (d), where # 7 is the polar 
angle of the radiative photon in the J/ip rest frame, and 9^ and (p> n are the polar angle and azimuthal angle of 77 
in the t]T] helicity frame. The PWA results provide a good description of data. To illustrate the contributions from 
each component, the PWA projections for each specific resonance are plotted (Fig. 3 (a)-(f): /o(1500), /o(1710), 
/ (2100), ./2( 1525 )> /2(1810), / 2 (2340)), 0++ phase space (Fig. 3 (g)), total 0++ component (Fig. 3 (h)) and total 
2 ++ component (Fig. 3 (i)), where the dots with error bars are data with the background events subtracted and the 
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TABLE I: Summary of the PWA results, including the masses and widths for resonances, branching ratios of J/tp — > 7X, as 
well as the significance. The first errors are statistical and the second ones are systematic. The statistic significances here are 
obtained according to the changes of the log likelihood. 



Resonance Mass(MeV/c 2 ) Width(MeV/c 2 ) B{ J/tp ->■ 7X -> 77777) Significance 



/o(1500) 


14681^ 


1 oc+41+28 
iJO -26-100 


(i-65iS:§;t°i:S) x 


-5 


8.2 


<T 


/o(mo) 


1759±6t2s 


172±10l?e 


(2.35±°;l?± ; 2 7 t) x 10" 


-4 


25.0 


(T 


/o(2100) 


2081±13l36 


OV o+27+70 
zl J -24-23 


(i.i3i° :?ro: 6 2 4 8) x 10- 


-4 


13.9 


(7 


/I (1525) 


1513±5±* 


-r + 12+16 
' d -10-8 


(3.42±°;S±i:S5) x 10- 


-5 


11.0 


rT 


/ 2 (1810) 


1822l£t? 7 


99Q+52+88 
■"''-42-155 


(5.40±g;~±i;S) x 10- 


-5 


6.4 


O 


/ 2 (2340) 


nq K 9+31+140 
^°"^-30-63 


oq/t+62+165 
°°*-54-100 


(5.60«;^ 2 2 ; 3 07 ) x 10" 


-5 


7.6 


(T 



solid histograms are the projections of the PWA for the specific components. 

1. Scalar components 

The histogram in Fig. 3 (h) shows the contribution of all the scalar components, where the dominant ones are 
from / (1710) and / (2100). For the / (1710) meson, the PWA gives a mass of 1759±6 MeV/c 2 and a width of 
172±10 McV/c 2 with a statistical significance of 25a; the mass and width are consistent with those obtained from 
J/tp -t jKK [23] and J/tp -> ^tttt [24] at BESII. The /o(2100) is observed with a statistical significance of 13.9c-, and 
its mass and width are determined to be 2081±13 MeV/c 2 and 273^ 2 4 MeV/c 2 , respectively, which are in agreement 
with previous measurements [25-28]. The product branching fractions of the /o(1710) and /o(2100) are measured to 
be B(J/tp -> 7/o(1710) -4 77777) = (2.35+g;^) x 10~ 4 and B(J/tp 7/ (2100) -> 7W )=(1.13±^;?§) x 10" 4 , where the 
errors are statistical only. 

The /o(1500) is observed with a statistical significance of 8.2a, but its production rate, B(J/tp — > 7/o(1500) — > 
77777) =(1.651q Jj) x 10 -5 , is about one order of magnitude lower than that of /o(1710) and /o(2100) since its dom- 
inant decay modes are 4-7T and 7T7t [16]. The mass and width obtained from the global fit are 1468^15 MeV/c 2 and 
136^26 MeV/c 2 , respectively, which are consistent with the BESII measurements in J/tp — > ^tttt [24]. 

The first experimental evidence for the / (1790) (M=1790t 4 ° MeV/c 2 and r=270ig° McV/c 2 ) was observed in 
J/tp — > 4>tttt [29]. Of interest is that no evidence was observed in J/tp — > <pKK [29]. In this analysis, if the dominate 
/o(1710) in the basic solution is replaced with /o(1790), the log likelihood is worsen by 30. If the /q(1790) is included 
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FIG. 2: Comparisons between data and PWA fit projections: (a) the invariant mass spectrum of r/rj, (b)-(c) the polar angle 
of the radiative photon in the J /if) rest frame and rj in the r\r\ helicity frame, and (d) the azimuthal angle of r\ in the fir? 
helicity frame. The black dots with error bars are data with background subtracted, and the solid histograms show the PWA 
projections. 

as an additional resonance in the fit, the significance of /o(1790) is only 1.8<7, which indicates that the not coupled 
/o(1790) is either suppressed in radiative decays or not coupled strongly to r/rj. 

To evaluate the contributions from other scalar mesons, ,/o(1370), /o(2020), /o(2200) and /o(2330), the PWA was 
performed including them, and none of them has significance greater than 5. Oct. Therefore, they are not included in 
the basic solution. 



2. Tensor components 

The total contribution from the tensor components is shown as the histogram in Fig. 3 (i), where the peak around 
1.5 GeV/c 2 is dominated by the well-established resonance / 2 (1525) and the tensor components contributing to the 
bump around 2.1 GeV/c 2 arc from / 2 (1810) and / 2 (2340). The fitted mass and width of / 2 (1525) arc 1513±5 McV/c 2 
and 751^ MeV/c 2 , respectively, which arc consistent with the world average values [16], and the product branching 
fraction is calculated to be B(J/ip — > 7/ 2 (1525) jrjr])=(3A2'to'^i) x 1CT 5 . If / 2 (1525) is replaced with another 
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FIG. 3: Contribution of the components, (a) /o(1500), (b) /o(1710), (c) /o(2100), (d) / 2 (1525), (e) / 2 (1810), (f) / 2 (2340), (g) 
++ phase space, (h) total ++ component, and (i) total 2 ++ component. The dots with error bars are data with background 
subtracted, and the solid histograms are the projection of the PWA. 

tensor meson close to 1.5 GeV/c 2 , /2(1565), the log likelihood is worse by 18. The PWA is also performed including 
/2(1565) as an additional resonance, and its statistical significance is only 2. Oct. 

The global fit shows that there is a tensor component around 1.8 GeV/c 2 with a statistical significance of 6.4ct, 
and its mass and width are determined to be 1822±24 MeV/c 2 and 229±H MeV/c 2 , respectively, which is likely to be 
the /2(1810). However the changes of the log likelihood value are only 0.8 or 0.7, if we replace it with the /2(1910) 
or /2(1950), respectively, using the world average values for their masses and widths [16], which indicates that we 
cannot distinguish it from /2(1810), /2(1910) and /2(1950) with the present statistics. In this analysis, this tensor 
component is denoted as /2(1810), and the ambiguous assignment of /2(1810) or /2(1950) is considered as a source 
of systematic error. 

To investigate contributions from other possible tensor resonances, /2(2010), /2(2150), /j(2220), /2(2300) and 
/2(2340), the fits were performed with alternative combinations, and the statistical significances of /2(2010), /2(2150) 
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and /j(2220) are all less than 5.0c, and the best fit favors the presence of /2(2340) (the statistical significance is 
7.6a) with a mass of 2362+^ MeV/c 2 , a width of 334^4 MeV/c 2 , and a product branching fraction of B(J/ip -> 
7/ 2 (2340) -> 7»m)=(5.60j:oj§) x 10 -5 . Since the mass of / 2 (2300) is close to / 2 (2340), an attempt was made to 
replace / 2 (2340) with / 2 (2300) by fixing its mass and width to those in PDG [16], and the log likelihood value is 
worse by 15. The narrow /j(2220) (also known as £(2230)), which was reported by Marklll [30] and BES [31], is also 
studied. In this analysis no evident narrow peak around 2.2 GeV/c 2 over the broad bump is observed in the rjrj mass 
spectrum shown in Fig. 1 (c). When the /,/(2220) is included in the PWA, the statistical significance is found to be 
0.4cr. 

3. Non-resonant contribution and J/ip — > <pr] 

In the analysis, the non-resonant contribution in the decay J/ip — > 77777 is described with ++ phase space, with 
a statistical significance of 12.4cr, and the product branching fraction is calculated to be (l-47ig o2) x 10~ 4 - An 
alternative fit is made by replacing the ++ phase space with 2 ++ phase space, and the log likelihood value is worsen 
by 30. In addition to ++ phase space, a 2 ++ phase space component to describe the non-resonant contribution was 
used, but the significance of the additional 2 ++ phase space is less than 4. Oct. The impact from the uncertainty of the 
non-resonant contribution is taken as a source of systematic error. 

V. SYSTEMATIC ERROR 

The systematic sources and their corresponding contributions to the measurements of mass, width and branching 
fractions are described below. 

• Background uncertainty. The background events estimated with the 77 mass sidebands are included in the global 
fit with negative weights. To estimate the systematic error, the global fit was done with background events from 
different 77 mass sideband regions, and the changes of results are assigned as the systematic errors. 

• Uncertainty from extra components. As mentioned above, possible extra ++ , 2 ++ and 4 ++ components with 
low significance were removed from the global fit. The changes of results caused by including them in the global 
fit are assigned as the systematic errors. 

• Uncertainty from resonance parameters. To estimate the impact of one specific resonance on the others, the 
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optimized mass and width of each resonance were varied by one standard deviation (statistical error only), 
and then the global fit was redone. The differences between the results with and without the variation of 
the resonance parameters are assigned as the systematic errors. As discussed above, the /2(1810) cannot be 
distinguished from ^(WIO) or /2(1950). Therefore the fits were redone fixing the mass and width to be the 
values of /2(1810) and /2(1950) in PDG [16], respectively, and the maximum changes of the results are regarded 
as the systematic uncertainties. 

• To estimate the uncertainty from J/ip —> <j)ij, an alternative fit was performed without the contribution of 
J/ip — > <prj, and the changes of results arc taken as the systematic errors. 

• Mass resolution. In the global fit, the mass resolution is not considered to simplify the analysis. In order to 
estimate its possible impact on the fitted resonance parameters, a test was made by smearing the line shape of 
each resonance from the global fit with the corresponding mass resolution obtained from MC simulation. The 
impact on a resonance with a width greater than 100 MeV/c 2 is less than 2%, which is negligible compared 
with uncertainties from other sources. For / 2 (1525), the width is smeared by 5 MeV/c 2 with respect to the 
PWA result 75^]q (stat.) MeV/c 2 . The difference is considered as a source of systematic error to the width 
measurement. 

• Phase space description. The uncertainty from the description of the non-resonant contribution is estimated 
from an alternative fit by including both ++ and 2 ++ phase space. 

• Breit-Wigner formula. The changes of the fit results caused by replacing the constant width Breit-Wigner 
with a kinematic dependent width Brcit-Wigncr [32] are taken as the uncertainties from different resonance 
paramcterizations. 

In addition to the above systematic sources, the systematic errors from the event selection criteria, trigger efficiency and 
the number of J/ip events, which are summarized in Table II, are also included in the branching fraction measurements. 

• Photon detection. For the decay mode analyzed in this paper, five photons are involved in the final states. The 
uncertainty due to photon detection and photon conversion is 1% per photon. This is determined from studies 
of photon detection efficiencies in well understood decays such as J/ip — > p°7r° and study of photon conversion 
via e + e~ — !> 77 [33, 34]. 
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TABLE II: Summary of the systematic errors from the event selection. 

Error sources Systematic error(%) 

Photon efficiency 5.0 

Kinematic fit 6.5 

r\ selection 0.8 

Number of J/ip events 1.24 

Total 8.3 



• Trigger efficiency. The trigger efficiency of the BESIII detector was found to be very close to 100% from studies 
using different samples selected from J/ip or ip{2S) decays. Therefore, the trigger efficiency is assumed to be 
100% in the calculation of the branching fractions, and the systematic error from this source is neglected. 

• Kinematic fit and rj selection. The systematic error from the kinematic fit is studied with the clean channel 
ip' —> 7XcO (Xco — > 1717)1 a s described in Ref. [34]. The efficiency is defined as the ratio of Xco yield with and 
without the kinematic requirement of xic ^ ^0, where the \c0 yield is obtained by fitting the r\r\ mass spectrum 
with the MC signal shape and a second-order polynomial. The difference between data and MC simulation, 
6.5%, is taken to be the systematic error. Similarly the systematic error from 77 selection criteria is estimated 
to be 0.8%. 

• Number of J/ip events. In the calculation of branching fractions, the number of J/ip events, (225.3±2.8) x 10 6 [9], 
determined from J ftp inclusive hadronic decays, was used, and its uncertainty, 1.24%, is taken as the systematic 
error. 

The systematic error sources and their contributions studied above are all summarized in Table III and Table IV, in 
which the systematic error from event selection includes the contributions from photon detection efficiency, kinematic 
fit, 77 selection and the number of J/tp events listed in Table II. The total systematic error is the sum of them added 
in quadrature. 

VI. SUMMARY 

Using 225 million J/ip events collected with the BESIII detector, a PWA of J/ip — > 77777 has been performed, and the 
results are summarized in Table I. The scalar contributions are mainly from /o(1500), /o(1710) and /o(2100), while 
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TABLE III: Summary of the systematic error sources and their corresponding contributions to masses and widths of the 
resonances X (MeV/c 2 ), which are denoted as AM and AF, respectively. 



Systematic error 


/o(lBOO) 


/a (1525) 


/o(1710) 


Mil 


510) 


/o(2100) 


/ 2 (2340) 




AM 


Ar 


AM 


Ar 


AM 


Ar 


AM 


Ar 


AM 


Ar 


AM 


Ar 


Background uncertainty 


+ 18 
-8 


-46 


+ 1 


+6 
— 4 


+ n 

— 13 


+ 19 
-8 


+ 55 
-39 


+ 61 
— 4 


+ 19 
-32 


+38 


+93 


+ 43 
— 41 


Extra resonances 


+9 
-73 


+ 7 
-84 


+ 2 
-9 


+ 10 
-4 


+8 

-18 


+ 19 
-11 


+ 19 
-36 


+ 16 
-141 


+ 7 
-9 


+ 50 
-4 


+ 76 
-9 


+ 7 
-2 


Resonance parameters 


+ 11 
-12 


+ 27 
-25 


+3 
—4 


+ 8 


+4 
— 11 


+ 17 
-7 


+24 
-18 


+ 61 
-26 


+ 13 
— 12 


+26 
-23 


+ 55 
-62 


+ 157 
-87 


J/tp -+ 4>TI 


-1 


-10 


(1 


+ 2 





+3 


+ 11 


+3 


-4 


+ 12 


+ 16 


+26 


Phase space description 





+3 


-1 


-1 


-1 





-11 


+7 


-3 


+ 10 


+38 


-21 


Brcit-Wigncr formula 





-6 


+ 1 


+6 


-5 


-5 


+ 15 


-58 


+2 


+4 


+20 


-18 


Total 


+ 23 
-74 


+ 28 
-100 


+4 
— 10 


+ 16 
— 8 


+ 14 
— 25 


+32 
— 16 


+66 
-57 


+ 88 
-155 


+ 24 
-36 


+ 70 
-23 


+ 140 
-63 


+ 165 
-100 



TABLE IV: Summary of the systematic error sources and their corresponding contributions to the branching fractions of 
J /ip — > jX, X — > rjt] (%), which are denoted as AB. 

Systematic error AS(/ (1500)) AS(/2(1525)) AS(/ (1710)) AS(/ 2 (1810)) AS(/ (2100)) AS(/ 2 (2340)) 



Event selection 


±8.3 


±8.3 


±8.3 


±8.3 


±8.3 


±8.3 


Background uncertainty 


+20.6 
-46.1 


+23.5 
-34.8 


+35.4 
-15.1 


+5.1 

-34.5 


+46.8 
-9.7 


+24.9 
-0.9 


Extra resonances 


+ 11.1 
-56.3 


+21.9 


+33.2 
-23.9 


+20.3 
-19.0 


+26.9 
-16.2 


+6.0 
-27.8 


Resonance parameters 


+ 18.0 

-41.6 


+21.6 
-12.0 


+ 17.0 
-8.0 


+58.2 
-14.4 


+ 11.3 
-11.1 


+24.3 
-19.3 


J/ip -> <Pv 


-7.1 


+0.6 


+7.6 


-6.4 


±8.7 


±5.0 


Phase space description 


-1.6 


-3.2 


-0.5 


±10.7 


-1.0 


+21.1 


Breit-Wigner formula 


-6.3 


+6.8 


-8.4 


-4.9 


-7.4 


-12.5 


Total 


+30.7 


+40.2 


+52.6 


+63.3 


+56.5 


+42.3 


-84.8 


-37.9 


-31.7 


-43.5 


-24.6 


-37.0 



no evident contributions from /o(1370) and /o(1790) are seen. Recently, the production rate of the pure gauge scalar 
glueball in J/tp radiative decays predicted by the lattice QCD [35] was found to be compatible with the production 
rate of J/tp radiative decays to /o(1710); this suggests that /o(1710) has a larger overlap with the glueball compared 
to other glueball candidates (eg. /o(1500)). In this analysis, the production rate of /o(1710) and /o(2100) are both 
about one order of magnitude larger than that of the /o(1500) and no clear evidence is found for /o(1370), which are 
both consistent with, at least not contrary to, lattice QCD predictions. 

Studies using data from pp annihilation [27, 28] show that the /o(2100) has strong coupling to 7777, but much weaker 
to tttt, which indicates an exotic /o(2100) decay pattern. Searching for more decay modes of /o(2100) in J/ip radiative 
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decays may help to clarify its nature. 

The tensor components, which are dominantly from / 2 (1525), /2(1810) and /2(2340), also have a large contribution 
in J ftp — > ^rj-q decays. The significant contribution from / 2 (1525) is shown as a clear peak in the 7777 mass spectrum; 
a tensor component exists in the mass region from 1.8 GeV/c 2 to 2 GeV/c 2 , although we cannot distinguish /2(1810) 
from /2(1910) or /2(1950); and the PWA requires a strong contribution from / 2 (2340), although the possibility of 
/2(2300) cannot be ruled out. For the narrow /j(2220), no evident peak is observed in the 7777 mass spectrum. We 
have also tried to add it in the analysis, but its statistical significance is quite small, just OAa. 
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